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One of the main open problems in the physics of A-hypernuclei is the lack of a sound theoretical 
interpretation of the large experimental values for the ratio T„/Vp = r(An — > nn)/r{Ap np). 
To approach the problem, we have incorporated a one-meson-exchange model for the A7V nN 
transition in finite nuclei in an intranuclear cascade code for the calculation of single and double- 
coincidence nucleon distributions corresponding to the non-mesonic weak decay of ^He and ^^C. 
Due to the elimination of interferences, two-nucleon coincidences are expected to give a cleaner 
determination of Vn/Vp than single-nucleon observables. Single-nucleon distributions are found 
to be less sensitive to variations of r^/Fp than double-coincidence spectra. The comparison of 
our results with preliminary KEK coincidence data allows us to extract a F„/Fp ratio for ^He of 
0.39 ± 0.11 when multinucleon induced channels are omitted. 



PACS numbers: 21.80.+a, 13.30.Eg, 13.75.Ev 



I. INTRODUCTION 



An old challenge of hypernuclear decay studies has 
been to secure the "elusive" theoretical explanation of 
the large experimental values of the ratio between the 
neutron- and proton-induced non-mesonic decay rates, 
P^/Pp = P(An nn)/T{Kp np) [1,2]. Indeed, the 
calculations underestimate the data for all considered 
hypernuclei. Moreover, in the experiments performed 
until now it has not been possible to distinguish be- 
tween nucleons produced by the one-body induced chan- 
nel, AA^ nN, and the two-body induced mechanism, 
ANN nNN, which is expected to be non-negligible 
and thus important for the determination of P^/Pp. 

Because of its strong tensor component, the one-pion- 
exchange (OPE) model with the A/ — 1/2 isospin 
rule supplies very small ratios, typically in the interval 
0.05 ~ 0.20. On the contrary, the OPE description can 
reproduce the total non-mesonic decay rates observed for 
light and medium hypernuclei. 

Other interaction mechanisms might then be necessary 
to correct for the overestimation of Pp and the underes- 
timation of P„ characteristic of the OPE. Those which 
have been studied extensively in the literature are the fol- 
lowing ones: i) the inclusion in the AA'^ —f nN transition 
potential of mesons heavier than the pion (also includ- 
ing the exchange of correlated or uncorrelated two-pions) 
[3-7] ; ii) the inclusion of interaction terms that explicitly 
violate the A/ =1/2 rule [1,8,9]; iii) the inclusion of the 
two-body induced decay mechanism [10-13] and iv) the 
description of the short range AA^ nN transition in 
terms of quark degrees of freedom [14,15], which auto- 
matically introduces A/ = 3/2 contributions. 

Recent progress has been made on the subject, 
(i) On the one hand, a few calculations [5-7,14] with 



AA^ — !■ nN transition potentials including heavy-meson— 
exchange and/or direct quark contributions obtained ra- 
tios more in agreement with data, without providing, 
nevertheless, an explanation of the origin of the puzzle 
[1]. In particular, these calculations found a reduction 
of the proton-induced decay width due to the opposite 
sign of the tensor component of if-exchange with respect 
to the one for 7r-exchange. Moreover, the parity violat- 
ing AA^('^S'i) — > nN{^Pi) transition, which contributes 
to both the n- and p-induced processes, is considerably 
enhanced by A'-exchange and direct quark mechanisms 
and tends to increase P„/Pp [6,14]. Very recently, the 
AA^ — > nN interaction has been studied within an ef- 
fective field theory framework [16]. The decay of s- 
and p-shell hypernuclei was approached following the 
same formalism as in Refs. [4,6], but the weak transi- 
tion was described in terms of OPE, one-kaon-exchange 
and [AS"! = 1 four-fermion contact terms. The results 
obtained in Ref. [16] are very encouraging and open a new 
door for systematic studies of hypernuclear weak decay 
based on effective field theory descriptions, 
(ii) On the other hand, an error in the computer program 
employed in Ref. [17] to evaluate the single-nucleon en- 
ergy spectra from non-mesonic decay has been detected 
and corrected in Ref. [18]. It consisted in the underes- 
timation, by a factor ten, of the nucleon-nucleon colli- 
sion probabilities entering the intranuclear cascade cal- 
culation. The correction of such an error leads to quite 
different spectral shapes and made it possible to repro- 
duce old experimental data for j^C [19,20] even with 
a vanishing value of Pn/Pp (which is a free parameter 
in the polarization propagator model of Refs. [17,18]). 
However, when compared with the recent proton energy 
spectra measured by KEK-E307 [21], the corrected dis- 
tributions still provide a quite large value of the ratio: 
P„/Pp = 0.87 ± 0.23 for ^^^c [22], which is incompati- 



1 



ble with pure theoretical estimations, ranging from 0.1 
to 0.5. 

In the light of these recent developments and of new 

experiments at KEK [23], FINUDA [24] and BNL [25], it 
is important to develop different theoretical approaches 
and strategies for the determination of Tn/Tp. In a pre- 
vious Letter [26] we discussed some results of an evalu- 
ation of nucleon-nucleon coincidence distribution in the 
non-mesonic weak decay of ^He and hypernuclei. 
The calculations were performed by combining a one- 
meson-exchange model describing one-nucleon induced 
weak decays in finite nuclei with an intranuclear cascade 
code taking into account the nuclcon final state interac- 
tions. The two-nucleon induced channel was also taken 
into account, treating the nuclear finite size effects by 
means of a local density approximation scheme. 

In the present paper we discuss the nucleon correla- 
tion observables in a more systematic way. In order to 
stress the importance of the coincidence analysis in con- 
nection with the determination of r^/Fp, we also discuss 
single-nucleon distributions. In principle, the correla- 
tion observables permit a cleaner extraction of F„/Fp 
from data than single-nucleon observables. This is due 
to the elimination of interference terms between n- and 
p-induced decays [1], which arc unavoidable in experi- 
mental data and cannot be taken into account by the 
Monte Carlo methods usually employed to simulate the 
nucleon propagation through the residual nucleus. We 
also perform a weak interaction model independent anal- 
ysis to extract an estimate for Fn/Fp using preliminary 
results from KEK [23,27] on two-nucleon angular and 
energy correlations. The resulting F„/Fp values for ^He 
turn out to be substantially smaller than those obtained 
in the past [20,28] from single coincidence analyses and 
fall within the predictions of recent theoretical studies 
[6,7,14]. 

The work is organized as follows. In Sec. II we give 
an outline of the models employed to describe the non- 
mesonic weak decay. In Sec. Ill we discuss the main 
features of the Monte Carlo simulation accounting for 
the nucleon propagation inside the residual nucleus. The 
purpose of Sec. IV is to discuss the reasons why multi- 
nucleon coincidence analyses should be preferred over 
single-nucleon studies in order to extract F„/Fp from 
data. Our results are discussed in Sec. V and the conclu- 
sions are given in Sec. VI. 



II. MODELS FOR THE WEAK DECAY 



A. One-nucleon induced decay: finite nucleus 
approach 

The one-nucleon induced non-mesonic decay rate can 
be written as: 
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where the initial hypernucleus, of mass Mh, is assumed 
to be at rest and the quantities En and Ei^2 are the 
total energy of the residual nuclear system and the to- 
tal asymptotic energies of the emitted nucleons, respec- 
tively. The integration variables Px = ki + k2 and 
kr = {ki — k2)/2 arc the total and relative momenta of 
the two outgoing nucleons. The momentum conserving 
delta function has been used to integrate out the mo- 
mentum of the residual nucleus, k^ ~ —Pt- The sum, 
together with the factor 1/(2J -|- 1), indicates an average 
over the initial hypernucleus total spin projections, Mj, 
and a sum over all quantum numbers of the residual sys- 
tem, {R}, as well as over the spin and isospin projections 
of the outgoing nucleons, {1} and {2}. Finally: 
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is the amplitude for the transition from an initial hyper- 
nuclear state into a final state which is factorized into 
an antisymmetrized two nuclcon state and a residual nu- 
clear state The two-nucleon state is characterized 
by the total momentum Pp, the relative momentum kr, 
the spin and spin projection S*, Ms and the isospin and 
isospin projection T, T3. Finally, Oan^uN is a two-body 
transition operator acting on all possible AA'' pairs. 

Working in a coupled two body spin and isospin basis, 
the non-mesonic decay rate can be written as: 

Ti = F„ -|- Fp, 

where F„ and Fp stand for the neutron (An nn) 
and proton-induced (Ap np) decay rate, respectively. 
They are given by {N = n,p): 
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where t^p = 1/2 and t^n = —1/2. The elementary ampli- 
tude tAN^nN accounts for the transition from an initial 
AA^ state with spin (isospin) 5*0 (To) to a final antisym- 
metric nN state with spin (isospin) S (T). It can be 
written in terms of other elementary amplitudes which 
depend on center of mass ("i?") and relative ("r") or- 
bital angular momentum quantum numbers of the AN 
and nN systems: 
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where the dependence on the spin and isospin quantum 
numbers has to be understood. In Eq. (2), the coeffi- 
cients X{NrLr, NrLr, IaIn) are the well known Moshin- 
sky brackets, while: 
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Here, Va.rir) stands for the one-meson-exchange weak 
potential, which depends on the relative distance between 
the interacting A and nucleon as well as on their spin and 
isospin quantum numbers. Moreover, $5v'l (^/(V26)) 
and $^^^^(.R/(6/\/2)) are the relative and center-of- 
mass harmonic oscillator wave functions describing the 
AA^ system, while (r) is the relative wave function 
of the nN final state. Further details regarding notation 
can be found in Ref. [4]. 

In order to study nucleon-nucleon coincidence spectra, 
it is convenient to introduce differential decay widths de- 
pending on the ccnter-of-mass coordinate R, the cosinus 
of the angle between nucleon 1 and 2, cos0i2, and the 
modulus of the momentum of one of the outgoing nucle- 
ons, say fci. For this purpose, taking into account Eqs. (2) 
and (3), let us rewrite the decay rate of Eq. (1) in the 
following schematic way: 
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where Ui and Cj include Clebsh-Gordan coefficients and 
other factors that depend on various quantum numbers, 
while Aij{R, Pt ,kr) denote AA^ nN weak decay am- 
plitudes. Changing variables from PT,kr to k\,k2, us- 
ing the energy conserving delta function and imposing 
rotational invariance, it is possible to substitute, in the 
integrand of Eq. (4), the kr and Pt dependence by a de- 
pendence on k\ and cos ^12. Without performing the R, 
ki, and cos0i2 integrations, one can then write n- and 
p-stimulated differential decay rates: 



rN{R, fci,cos6ii2) = 
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appropriate for obtaining the distributions of the weak 
decay nuclcons required as input of the intranuclear cas- 
cade calculation. As explained in the following Section, 
the intranuclear cascade code allows then the primary 
nucleons to change energy, direction and charge, exciting 
other secondary nucleons that are also followed as they 
move through the nucleus. 

Note that by considering the primary nucleons as 
emerging from the same point R in space, we have im- 
plicitly assumed that the weak transition proceeds as 
a point interaction. This is implied by the shape of 
the regularized potentials of the model of Refs. [4,6] — 
shown in Ref. [29] — , which peak at relative distances 
r =\ri — < 0.5 fm. 

Finally, we recall that in our calculations the one- 
meson-exchange (OME) weak transition potential enter- 
ing Eq. (3) contains the exchange of p, K, K* , lu and 
r] mesons in addition to the pion [6]. The strong cou- 
plings and strong final state interactions acting between 
the weak decay nucleons are taken into account by us- 
ing a scattering nN wave function from the Lippmann- 
Schwinger (T-matrix) equation obtained with NSC97 
(versions "a" and "f " ) potentials [30] . The corresponding 
decay rates for ^He and ^^C are listed in Table I (OMEa 
and OMEf) together with the one-pion-exchange predic- 
tions (OPE). 



B. Two— nucleon induced decay: local density 
approximation 

The differential and total decay rates for the two- 
nucleon induced process, Anp nnp, are calculated with 
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the polarization propagator method in local density ap- 
proximation (LDA) ofRefs. [11,12]. In such a calculation, 
the simple OPE mechanism, supplemented by strong 
AiV and nN short range correlations (given in terms of 
phenomenological Landau functions) described the weak 
transition process. The two-nucleon induced differential 
decay width r2(-R, fci, cos 6*12) is obtained from the two 
particle two-hole (2p2h) contributions to the pion self- 
energy. As explained in Ref. [11], these contributions are 
derived in a phenomenological way from a fit to pionic 
atoms, conveniently modified by the Lorentz-Lorenz cor- 
rection and extended to all kinematical regions using the 
appropriate phase-space. The intranuclear cascade code 
then follows the fate of the two nucleons excited by these 
2p2h components, as well as that of the third nucleon 
emitted at the AA'^tt vertex, as they move through the 
nucleus. 

In the present calculation, the distributions of the weak 
decay nucleons and the value of r2 have been properly 
scaled to maintain the ratio unchanged: we then 

use Fa/ri = (ra/ri)^^^ = 0.20 for ^He and 0.25 for 
\^C. These values of {T2/Ti)^^^ have been obtained with 
the latest update of parameters given in Ref. [12] after 
correcting a small (conceptual) error in the implemen- 
tation of data on the P- wave pion nucleus optical po- 
tential. For the hypernuclei treated in the present paper, 
I^He and j^^C, such a correction slightly decreases the val- 
ues of Fi while increasing F2 by about 20%. 



III. INTRANUCLEAR CASCADE SIMULATION 

In their way out of the nucleus, the weak decay (i.e., 
primary) nucleons continuously change energy, direction 
and charge due to collisions with other nucleons. As a 
consequence, secondary nucleons arc also emitted. 

We simulate the nucleon propagation inside the resid- 
ual nucleus with the Monte Carlo code ofRefs. [17,18]. A 
random number generator determines the decay channel, 
n-, p- or two-nucleon induced, according to the ratios 
F„/Fp and F2/F1 predicted by our finite nucleus and 
LDA approaches. Positions, momenta and charges of the 
weak decay nucleons are selected by the same random 
number generator, according to the corresponding prob- 
ability distributions given by the finite nucleus and LDA 
calculations. 

For neutron-, proton- and two-nucleon induced de- 
cays, the discussion of Section II allows us to in- 
troduce the differential decay rates F„(.R, /ci, cos0i2), 
Tp{R, fci, cos^i2) and F2(^, fci, cos^i2) supplying the n-, 
p- and two- nucleon stimulated total rates F„, Fp and F2 
through the following relation: 

Ti = I (fki / dcos6li2 / (fRTi{R,ki,cosei2)- 



After they are produced as explained, the pri- 
mary nucleons move under a local potential Vn{R) = 
-k%iR)/2mN, where kp^iR) = [Sn^pNiR)]^/^ (N = 
n,p) is the local nucleon Fermi momentum correspond- 
ing to the nucleon density pn{R)- The primary nucleons 
also collide with other nucleons of the medium according 
to free space nucleon nucleon cross sections [31] properly 
corrected to take into account the Pauli blocking effect. 
For further details concerning the intranuclear cascade 
calculation see Ref. [17]. Each Monte Carlo event will 
then end with a certain number of nucleons which leave 
the nucleus along defined directions and with defined en- 
ergies. One can then select the outgoing nucleons and 
store them in different ways, as we shall do in the discus- 
sion of Section V. 

We are aware of the fact that accounting for nucleon 
final state interactions effects in light residual nuclei (as 
those required to treat a He) through Monte Carlo tech- 
niques is questionable. However, realistic calculations 
for few-body scattering states have been performed up 
to 3 nucleons only (see for instance Ref. [32]). For the 
hypernuclear non-mesonic decay problem, only the case 
of has been treated exactly [29,33]. Although one 
might attempt three-cluster type calculations, this goes 
beyond the scope of the present work. For this reason, 
the results we present for ^He should be considered less 
realistic than the corresponding ones for j^^C. 

IV. SINGLE-NUCLEON VS 
NUCLEON-COINCIDENCE ANALYSES 

In this Section we discuss the reason why multi- 
nucleon coincidence studies should be preferred over 
single nucleon analyses when the purpose is the deter- 
mination of r„/Fp . The argument for this explanation 
is based on the elimination of quantum mechanical inter- 
ferences between n and p-stimulated weak decays [1]. 

Let us first note that the nucleons originating from 
n- and p induced processes are added incoherently (i.e., 
classically) in our intranuclear cascade calculation (for 
the moment we are then neglecting an analogous effect 
due to the two-nucleon stimulated decay channel). How- 
ever, for particular kinematics of the detected nucleons 
(for instance at low kinetic energies), an in principle 
possible quantum-mechanical interference effect between 
n- and ^(-induced channels should inevitably affect the 
observed distributions. Therefore, extracting the ratio 
F„/Fp from experimental data with the help of a classi- 
cal intranuclear cascade calculation may not be a clean 
task. 

To clarify better the issue, let us consider for instance 
an experiment (such as the majority of the experiments 
performed up to now) measuring single-proton kinetic 
energy spectra. The relevant quantity is then the num- 
ber of outgoing protons observed as a function of the 
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kinetic energy Tp. Schematically, this observable can be 

written as: 



where the various arguments have been suppressed and 
the number of nucleons; 



Np{Tp)(X (p(rp)|dFSlOwD|*H) 



a (p(Tp)|6FSi|nn, ^r) + (3 {p{Tp)\6vsi\np, 



(6) 



where \p{Tp)) represents a many-nucleon final state with 
a proton whose kinetic energy is Tp. Moreover, in Eq. (6) 

the action of the weak decay operator Owd = 0\n—ynn + 
OAp^np produced the superposition: 

Owd|*h) = Q!|nn,*ij) +/3|np,*ij/). 

Here |nn, ^'_r) ^j?/)) is a state with a nn (np) pri- 
mary pair moving inside a residual nucleus "^r {^w). 
Note that in the present schematic picture: r„ a |ap 
and Tp oc Since both transition amplitudes en- 

tering the last equality of Eq. (6) are in general non- 
vanishing, interference terms between n- and p-induced 
decays are expected to contribute to Np{Tp). An am- 
plitude {p{Tp)\OFSi\nn,'i^R) different from zero means 
that, due to nucleon final state interactions (FSI), a sec- 
ondary proton has a non vanishing probability to emerge 
from the nucleus with kinetic energy Tp even if the weak 
process was n-induced (i.e., without primary protons). 
While for high kinetic energies (say for Tp > 80 MeV) 
this amplitude is expected to be almost vanishing, as 
long as Tp decreases its contribution could produce an 
important interference effect (see the results discussed in 
Section VA). 

An interference-free observation would imply the mea- 
surement of all the quantum numbers of the final nucle- 
ons and residual nucleus. While this is an impossible 
experiment, what is certain is that the magnitude of the 
interference can be reduced if one measures in a more ac- 
curate way the final state. For this reason, two-nucleon 
coincidence observables are expected to be less affected 
by interferences than single-nucleon ones and thus more 
reliable for determining F^/Fp. If we consider the detec- 
tion of the kinetic energy correlation of a np pair emitted 
nearly back-to-back (say with cos^„p < —0.9): 



Ar„p(T„ + Tp,cosenp < -0.9) 



(7) 



cx 



a{n{T„),p{Tp), COS Onp < -0.9|6FSi|nn, -J-ij) 



+f3{n{Tn),p{Tp), cos 0np < -0.9|dFSi|np, , 

one expects an interference effect smaller than the one 
appearing in single-nucleon observations, i.e., in Eq. (6), 
especially when particular energy cuts are considered. 

To see this more explicitly, let us rewrite Eqs. (6) and 
(7) in the following way: 



Np = [ATA"-"" + ATAp-np] [1 + COS 0p] , (8) 

Nnp = " + A^^^-^"^] [1 + /„p cos cp^p] , (9) 



An—^nn 
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N^P^^P = ^ l/^l' \{np\dFsi\np,^R') 

have obvious meaning. Moreover: 
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and 



±An^nn 



(l>np~^^^) is the phase difference between the ampH- 



"^np 

tudes A^"^"" = y/A'A"- 



and AAp^"P 



AjAn—^nn . 
np 



and 



AAp^np 
np 



Note that the dis- 
tributions of Eq. (8) [Eq. (9)] are not affected by in- 
terference only when Ip cos = [/„p cos (jy^p = 0] . 
Since, as explained, we expect N^^^''" /N^p^^'p < 
j^An^nnjjqAp^np ^ ^ ^.j^gj^ ^jj^ smaller than Ip. 

The numerical results discussed in Sections V A and V B 
confirm this expectation. 

The same reasoning must be applied — and the pre- 
vious equations changed accordingly — once the two- 
nucleon induced decay mechanism is taken into account. 
Note, however, that our discussion in this Section have 
neglected the effect of the phase differences (j)p and (jjnp, 
which cannot be evaluated theoretically. A quite indirect 
estimate of these phases (and then of the real interfer- 
ence effects) can be obtained through the comparison of 
our calculated spectra (in which interference is not taken 
into account) with the experimental ones (which can be 
affected by interferences). 



V. RESULTS 



A. Single-nucleon spectra 



In Figs. 1 and 2 we show the proton kinetic energy 
spectra corresponding to the decay of ^He and ^^^C, re- 
spectively. Note the particular normalizations of the 
curves presented in this and in the following figures. 
The dashed curves correspond to the distributions of 
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the one- nuclcon induced primary protons: they have 
been obtained through intranuclear cascade calculations 
in which the one-nucleon stimulated weak decay nucle- 
ons move under the effect of the nuclear local potential, 
Vjv(-R) = —kp^{R)/2mN, without colliding with the nu- 
cleons of the medium. The inclusion of nucleon FSI pro- 
vide the result given by the dotted lines. The continuous 
lines correspond to the full calculation, i.e., once the two- 
nucleon induced channel is also included. The calcula- 
tions have been performed with the OMEf model, which 
predicts Tn/Tp = 0.46 for ^He and Tn/Tp = 0.34 for 
^C. The model OMEa supplies similar results both for 
the proton and neutron spectra. While the primary pro- 
ton distributions are very similar for the two hypernuclei 
(the one for j^^C is slightly broader), the full calculations 
clearly show that, as expected, FSI have a bigger effect 
for the heavier system. For ^•^C, FSI are so important 
that they completely smear out the maxima correspond- 
ing to the primary protons. Our predictions for the pro- 
ton spectra from j^C should be compared with the one 
measured at KEK-E307 [21]. Unfortunately, this is not 
possible since these data have not been corrected for the 
detector geometry and the energy losses occurring in the 
target and detector materials. 



0.4 p 
0.35 

0.3 
0.25 



Np 0.2 



0.15 

0.1 - 
0.05 
— 



IN without FSI 

IN+FSI 

1N+2N+FSI 



OMEf 



20 40 




100 120 140 



T (MeV) 



FIG. 2. Same as in Fig. 1 for i^C. 



Let us now introduce the number of nucleons of the 
type N (N = n,p) produced in n-induced (iV^^"), 
;)-induced (Nj^^) and two-nucleon induced (N^) de- 
cays. If we normalize these quantities per n-induced, 
p-induced, and two-nucleon-induced decay, respectively, 
the total number of nucleons of the type N is given by: 
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T (MeV) 

FIG. 1. Single-proton kinetic energy spectra for the 
non-mesonic decay of aH^- The dashed and dotted lines are 

normalized per one-nucleon induced decay (Fi = r„ -I- Fp), 
while the continuous line is normalized per non-mesonic de- 
cay (Fnm =Fi +F2). 



Nn 



ATlBn p 



(10) 



TyAn- 



-I- JVjy , 



Ap- 
N 



and Ar;^"f-""f have obvious 



where iV^"" 

meaning. All these nuclcon numbers can be consid- 
ered either as being functions of the nucleon kinetic 
energy, Nn{Tn), as it is the case of Figs. 1 and 2, 
or as the corresponding integrated quantities, Nn = 
JdTNNNiTN). By construction, NjP'', iVjfP and A^^^^ 
(AA"-"«, 7V^P^"f and 7V^"P-"»P) do not depend {do 
depend) on the interaction model employed to describe 
the weak decay. 



In Figs. 3 and 4, Nn, 

■Anp — >nnp 
N 



ArAn- 



N 



Ap^np 
N 



and 



N 

energy in the case of ^Hc. 



are shown as functions of the nucleon kinetic 
In Figs. 5 and 6 the same 
quantities are given for ^-^C. From the Nn and Np dis- 
tributions for ^Hc wc note that the maxima occurring 
at Tn — 75 MeV — mainly due to the kinematics of the 
weak decay nucleons (see the dashed line in Fig. 1 for 
protons) — are more pronounced for neutrons than for 
protons. Note that for any value of the ratio F^/Fp, the 
number of primary neutrons is indeed larger than the 
number of primary protons [see Eq. (11)]. Consequently, 
due to neutron-proton reactions, the proportion of sec- 
ondary protons in the proton spectrum is larger than the 
proportion of secondary neutrons in the neutron spec- 
trum. From Figs. 3 and 5 we also note that the fractions 
of protons emitted in netiiron-induced processes are quite 



6 



small. For values of Tp in the 30 40 MeV bin (i.e., just 
above the experimental threshold), 7V^«^»"/ArAp^"P ~ 
0.15 both for ^He and j^^C, which corresponds to an inter- 
ference term /p cos (^p = 0.67 cos (/)p in Eq. (8). Therefore, 
even if the number of protons from n-induced decays is 
expected to be rather small, the existence of such events 
could nevertheless produce a quite big interference ef- 
fect: unfortunately, the matter depends crucially on the 
phase which we cannot estimate theoretically. Prom 
Figs. 3 and 5 one also note that the two nuclcon induced 
mechanism could be responsible for an interference ef- 
fect in the proton energy spectra as important as the one 
which could be due to the n-induced channel. For single- 
neutron spectra these interference effects could be even 
larger. Indeed, from Figs. 3-6 it is clear that, as expected, 

jyAn^nn^j^Ap^np ^ jyAn^nn ^ j^Ap^np ^ ^ QnCC again 

this conclusion neglects the effect of the phases of the 
interfering amplitudes. 
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FIG. 4. Single-neutron kinetic energy spectra for the 
non-mesonic weak decay of aH^- The total spectrum N„ 
(normalized per non-mesonic weak decay) has been decom- 
posed in its components iV^"^"", N^p^"^ and N^^p^^^p 
according to Eq. (10). 
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FIG. 3. Single-proton kinetic energy spectra for the 
non-mesonic weak decay of aH©- The total spectrum Np 

(normalized per non-mesonic weak decay) has been decom- 
posed in its components N^"^"", N^p^'^p and N^rip^^np 
according to Eq. (10). 
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FIG. 5. Same as in Fig. 3 for i^C. 



The single-neutron spectrum for ^■^C observed in the 
KEK-E369 experiment [34] is well reproduced by our cal- 
culations. This is evident from Fig. 6, where we show 
results based on two models (OPE and OMEf) which 
predict quite different F^/Fp ratios. Unfortunately, the 
dependence of the neutron spectra on variations of F^/Fp 
is very weak (the same is true also for the proton spectra) 
and a precise extraction of the ratio from the KEK-E369 
distribution is not possible. We checked that an analo- 
gous calculation performed with the polarization propa- 
gator method in local density approximation of Ref. [18] 
supplies neutron spectra which reproduce the data of 
Fig. 6 with a per D.O.F. smaller than 1 when F^/Fp 
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(a free parameter in such kind of calculations) is chosen 
to lie in the range 0-^1.5 and data above r„ = 30 MeV are 
considered. The little sensitivity of 7V„ and Np to Tn/Tp 
is mainly due to the fact that these numbers arc normal- 
ized per non-mesonic weak decay. As a consequence, one 
should separately compare the complementary observ- 
able, r„ + Tp, with the experiment. For ]^G, our calcula- 
tions supply Tnm = Tn + rp + T2 = 1.25 {Tn + Tp) = 0.91 
or 0.69 when model OMEa or OMEf is used. These val- 
ues agree quite well with the KEK data 0.83 ± 0.11 of 
Ref. [35]. 



0.5 - total 

total OPE 

from n-induced 

4 - from p-induced 

from 2N-induced 

• KEK-E369 
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T (MeV) 

FIG. 6. Same as in Fig. 4 for fC. The OPE result is also 
shown. Experimental data are taken from Ref. [34] . 

The problem of the small sensitivity of Nn and Np to 
variations of Tn/Tp can be overcome if one concentrates 
on another single-nucleon observable. The ratio Tn/Tp is 
defined in terms of the ratio between the number of pri- 
mary weak decay neutrons and protons, N^^ and Np^, 
in the following way: 

Tn^l(N^_ \ 

Tp - 2 \N^'i V • 

Due to two-body induced decays and (especially) nucleon 
FSI, one expects the inequality: 

Y'p^liwp-')^^^^'''^-^^^''^^^ 

to be valid in a situation, such as the experimental one, 
in which particular intervals of variability of the neutron 
and proton kinetic energy, Ar„ and ATp, are employed 
in the determination of Nn and Np. 

As one can deduce from previous figures, not only Af„ 
and Np but also the ratio Nn/Np depends on AT„ and 
ATp. This is more evident from Table II, in which the 
function Ri defined by Eq. (12) is given for ^He and 
j^C, for different nucleon energy thresholds Tj^ and for 



the OPE, OMEa and OMEf models. For a given en- 
ergy threshold, Ri is closer to Tn/Tp for ^He than for 
j^C since FSI are larger in carbon. The ratio Nn/Np (or 
equivalently i?i) is more sensitive to variations of Tn/Tp 
(see the differences between the OPE, OMEa and OMEf 
calculations of Table II) than and Np separately. 
Moreover, Nn/Np is less affected by FSI than Nn and 
Np. Therefore, measurements of Nn/Np should permit 
to determine Tn/Tp with better precision. The recent 
KEK-E462 experiment has measured the ratio Nn/Np 
for %B.e: a preliminary analysis of the data supplies a 
value of i?i around 0.6 with a relative error of about 20% 
using nucleon energy thresholds of 50 and 60 MeV [36]. 
Our result of Table II corresponding to the OMEf cal- 
culation for Tjy = 60 MeV agree with this experimental 
determination. According to this comparison, the ratio 
Tn/Tp for ^He should be around the value (0.46) pre- 
dicted by the OMEf model. 

B. Double— coincidence spectra 

In Fig. 7 we report the opening angle distribution of np 
pairs emitted in the non-mesonic decay of ^He. Note the 

particular normalizations of the curves presented in this 
and in the following figures. Neglecting nucleon FSI and 
the two-nucleon induced channel, the angular correlation 
is strongly peaked at 9np = 180° and only 1 pair out of 
100 is emitted with an opening angle smaller than 115° 
(i.e., with cos 9np > -0.4). The effect of FSI is to de- 
crease the weak decay nucleon back-to-back distribution 
(i.e., for cos 9np < —0.9) of about 25% and populate, 
strongly, the spectrum for cos 6np > —0.8. The effect 
of the two-nucleon induced channel is moderate: practi- 
cally, it only increases the distribution (by about 20%) 
in the region with cos 0np > 0.4. When a nucleon kinetic 
energy cut of 30 MeV is applied, large part of the distri- 
bution is removed. This is due to FSI, which leads to a 
large amount of pairs (~ 70% of the total) in which at 
least one nucleon has kinetic energy below 30 MeV. 

In Fig. 8 we show the kinetic energy correlation of np 
pairs emitted in the decay of %ile. Neglecting nucleon 
FSI and the two-nucleon induced channel, the energy 
correlation is strongly peaked, as expected, at T„ + Tp ~ 
155 MeV. Indeed, the Q-value corresponding to the 
proton-induced three-body process ^He ^ "^H -I- n is 
153 MeV. The effect of the FSI is to decrease the back-to- 
back maximum and to populate, strongly, the spectrum 
for Tn+Tp < 140 MeV. The effect of the two-nucleon 
induced channel is only visible when T„ -|- Tp is below 
110 MeV, where it enhances the distribution. Once a 
30 MeV kinetic energy cut is applied, the distribution at 
small T„-|-Tp is considerably reduced for the same reason 
explained in the previous paragraph. 

The opening angle and kinetic energy correlations for 
nn pairs have essentially the same structure of the np 
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distributions showed in Figs. 7 and 8. For a discussion 
of the pp distributions and the different effect of FSI in 
\E.e and j^^C we refer to our previous paper [26]. 
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FIG. 7. Opening angle distribution of np pairs emitted 
in the non-mesonic decay of aHc. The continuous and 
dotted hues are normalized per one-nucleon induced decay 
(Fi = F„ + Fp), while the dashed and dot-dashed curves are 
normalized per non-mesonic decay (Fnm = Fi + F2). 



when Nnn and Nnp arc determined by employing particu- 
lar intervals of variability of the pair opening angle, A9i2, 
and the nucleon kinetic energies, AT„ and ATp. Actu- 
ally, as one can deduce from the figures wc discussed 
in Ref. [26], not only Nnn and Nnp but also the ratio 
Nnn/Nnp depends on A^i2, Ar„ and ATp. The discus- 
sion of Rcf. [26] also proves that the ratio Nnn /Nnp is 
much less sensitive to FSI effects and variations of the 
energy cuts and angular restrictions than Nnn and Nnp 
separately. 

The numbers of nucleon pairs Nj^n discussed up to 
this point are related to the corresponding quantities for 
the one nucleon {Nj^j^) and two- nucleon (A^^^) induced 
processes [the former (latter) being normalized per one- 
body (two-body) stimulated non-mesonic weak decay] 
via the following equation: 



N 



NN 



IB 

NN 



Fi + F2 



(14) 



where: 



N 



N 



IB 
NN 



Ara- 
NN 



-r ^^NN 



-1- ^^NN 



■nnp 



ArlBn p I ArlBp p 
-''JVJV '- n ^ -"'JViV P 



(15) 
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FIG. 8. Kinetic energy sum distribution of np pairs emit- 
ted in the non-mesonic decay of aHs. The continuous and 

dotted Unes are normalized per one-nucleon induced decay 
(Fi = F„ + Fp), while the dashed and dot-dashed curves are 
normalized per non-mesonic decay (Fnm = Fi -I-F2). 



The ratio F^/Fp is defined as the ratio between the 
number of weak decay nn and np pairs, Nnn and Nnp- 
However, due to two body induced decays and (espe- 
cially) nucleon FSI effects, one has: 



Tn 

F„ 



N 



wd 



N, 



Afwd 
np 



^ = i?2 [A^i2, AT„, ATp, T2] (13) 



and the remaining iV's have obvious meaning. Therefore, 
the quantities Nj^^, AT^^^ and Ar2B^ (iV^V""- -^^y"" 



and iv;)';;; 



Anp-tnnp-^ do not depend {do depend) on the inter- 
action model employed to describe the weak decay. 

In Table III we report the OPE, OMEa and OMEf pre- 
dictions for Nnp and its components iV^p"^"", N^p^"p 
and N^^P^'^'-'P in the case of ^^^^^ ^wo different open- 
ing angle regions and an energy threshold of 30 MeV 
have been considered. The contribution to Nnp of the 
n-induced decay channel is always smaller than or com- 
parable with the one coming from the two-nucleon in- 



duced mechanism: N, 



AT An — >nn /AjAp 
np I np 



An- 
np 



< 



J\jAnp- 
np 



'Unp 



The ratio 



^P is always smaller than 0.10, which cor- 
responds to Inp < 0.57 in Eq. (9). This shows that 
interferences in coincidence obscrvablos arc potentially 
smaller than in single-nucleon spectra but also that they 
might be non-negligible. Due to the less pronounced 
effects of FSI in ^Hc, smaller values of Inp have been 
obtained in this second case. 

In Fig. 9 we show the np pair opening angle distribu- 
tion in the case of PC. The total spectrum Nnp has been 



decomposed into the components Nnp^' , ^-np 



AfAnp 
np 



np 



N^P^^'P and 

'""'P. A nucleon energy threshold of 30 MeV has 
been used in the calculation. Fig. 10 corresponds to the 
kinetic energy correlation of np pairs: it is again for 
and = 30 MeV, but now only back-to-back angles 
(cos 9np < —0.8) have been taken into account. We note 
how both the n-induced and the two-nucleon induced de- 
cay processes give very small contributions to the total 
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distributions in Figs. 9 and 10. Nevertheless, these de- 
cay processes could produce non-negligible interference 
terms. To minimize this effect, one could consider, for 
instance, not only back-to-back angles but also nucleon 
kinetic energies in the interval 150-^ 170 MeV, for which 
we predict /„p = 0.18. Again, smaller values of /„p have 
been obtained for ^^He. 
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FIG. 9. Opening angle distributions of np pairs emitted in 
the non-mesonic decay of a^C. The total spectrum Nnp (nor- 
malized per non-mesonic weak decay) has been decomposed 
in its components iV^""*"", iV^^^"" and iV^"''^""^ accord- 
ing to Eqs. (14) and (15). 



old of 30 MeV. The results of the OMEa and OMEf 
models are in reasonable agreement with the preliminary 
KEK-E462 data for ^He [27]. 

This data {Nnn/Nnp = 0.44±0.11), which corresponds 
to an energy threshold of 30 MeV and cos Onp < —0.8, 
can be fitted using the 6 weak interaction model inde- 
pendent quantities AT^^", N^^, N^, Nl^^, AT^^p and 
N"^ entering Eqs. (14) and (15) and quoted in Table V. 
This can be achieved through the following relation: 



j^lBn , JY2Bi£ \ £« , ^IBp 
1 1 / 1 p 



.A^2bE2 
nn 

t 1 



^ihin , Jv2Bi2_ \ l_a ^IBp , ^2B 

np ' np y r ^ np ' np 



r„ 



(16) 



using Fn/Fp and F2/F1 as fitting parameters. 

In Fig. 11 (Fig. 12) we report the dependence of 
Nnn/Nnp foi ;\Hc (^^C) On the ratio Vn/Tp for four dif- 
ferent values of r2/Fi. Both figures correspond to the 
case with a nucleon energy threshold of 30 MeV and the 
angular restriction cos Onp < —0.8. For a given value of 
r2/Fi, Figs. 11 and 12 permit an immediate determina- 
tion of F^/Fp by a direct comparison with data for the 
observable Nnn/Nnp. 

Using the \YIq data Nnn/Nnp = 0.44 ± 0.11 and as- 
suming F2 = 0, from Eq. (16) we obtain the following 
fitted ratio: 



0.1 



N.. 



0.05 



total 

from n-induced 
from p-induced 
from 2N-induced 



OMEf 



T,T > 30 MeV 

n' p 

COS 0„p < -0.8 



^ (iHe) =0.39 ±0.11. 

Tp 



(17) 



By employing the value F2/F1 = 0.2 (i.e., the one ob- 
tained with the model of Ref. [12] and used in our calcu- 
lations), a 34% reduction of the ratio is predicted: 



(iHe) =0.26 ±0.11. 

-L T) 



(18) 
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FIG. 10. Kinetic energy correlations of np pairs emitted in 
the non-mesonic decay of a^C at cos 9np < —0.8. The total 
spectrum Nnp (normalized per non-mesonic weak decay) has 
been decomposed in its components N^p^"", N^p^"'' and 
j^Anp^nnp according to Eqs. (14) and (15). 



In Table IV the ratio Nnn/Nnp predicted by the OPE, 
OMEa and OMEf models for ^He and is given for two 
opening angle intervals and for a nucleon energy thresh- 



These values of F„/Fp are rather small if compared 
with previous determinations [20] (0.93 ± 0.55) [28] 
(1.97 ± 0.67) from single— nucleon energy spectra anal- 
yses. On the contrary, the results of Eqs. (17) and (18) 
are in agreement with the pure theoretical predictions of 
Refs. [6,7,14]. 

Forthcoming coincidence data from KEK and FIN- 
UDA could be directly compared with the results dis- 
cussed in this paper and will permit to achieve new de- 
terminations of the F„/Fp ratio and to establish the first 
constraints on r2/ri. 
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FIG. 11. Dependence of the ratio Nnn/Nnp on Fn/Fp 
and F2/F1 for iHe. The results correspond to a nucleon 
energy threshold of 30 McV and the angular restriction 
cos 9nn < —0.8. The horizontal lines show the preliminary 
KEK-E462 data of Ref. [27]. 
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FIG. 12. Same as in Fig. 11 for fC. 



VI. CONCLUSIONS 

In this work we have presented a calculation of sin- 
gle and double-coincidence observables for the nucleons 
emitted in the non-mesonic weak decay of A-hypemuclei. 
This has been possible by supplementing our OME weak 
interaction models with FSI through an intranuclear cas- 
cade calculation. 

The results of our calculations are in reasonable agree- 
ment with preliminary KEK-E462 data for the observ- 
able ratio Nnn/Nnp in 

We also pc^rform a weak interaction model independent 
analysis in which Tn/Tp and T2/Tx are considered as free 
parameters: the KEK |He data Nnn/Nnp = 0.44 ± 0.11 



is reproduced if r„/rp = 0.39 ± 0.11 and r2 = or 
r„/rp = 0.26 ± 0.11 and Fs/ri = 0.2. 

The extension of the present study to triple-nucleon 
coincidence is of interest both for the determination of 
r„/rp and to disentangle the effects of one- and two- 
nucleon induced decay channels. 

The values of T„/Tp wo obtain by fitting KEK coin- 
cidence data for ^He are in agreement with other re- 
cent theoretical evaluations. However, they are rather 
small if compared with the results of previous analy- 
ses from single-nucleon energy spectra. Actually, all the 
previous experimental analyses of single-nucleon spectra 
[19-22,28,37], supplemented in some cases by intranu- 
clear cascade calculations, derived r„/rp values in dis- 
agreement with pure theoretical predictions. In our opin- 
ion, the fact that our calculations reproduce coincidence 
data for values of Tn/Tp as small as 0.3 -i- 0.4 could sig- 
nal the existence of non-negligible interference effects be- 
tween the n- and p-induced channels in those old single- 
nucleon data. 

In conclusion, although further (theoretical and exper- 
imental) confirmation is needed, we think that our in- 
vestigation proves how the study of nucleon coincidence 
observables can offer a promising possibility to solve the 
longstanding puzzle on the Tn/Tp ratio. 
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TABLE I. Weak decay rates (in units of the free A decay 
width) predicted by Ref. [6] for ^He and ]^C. 
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TABLE II. Predictions for the quantity Ri of Eq. (12) for 
^Hc and corrcsi)on(liiig to difforont nucleoli thresholds 
and to the OPE, OMEa and OMEf hkkIoIs. 









(McV) 













30 


60 






OPE 


0.04 


0.13 


0.16 


0.09 


iHe 


OMEa 


0.15 


0.32 


0.39 


0.34 




OMEf 


0.19 


0.40 


0.49 


0.46 




OPE 


-0.06 


-0.01 


0.05 


0.08 


12^1 
A 


OMEa 


-0.02 


0.07 


0.19 


0.29 




OMEf 


-0.01 


0.09 


0.21 


0.34 



TABLE HI. OPE, OMEa and OMEf predictions for N^p 

and its components N^^^"" , N^^^"^ and iV^""^""'' (in- 
tegrated over all angles and for energies Tjv > 30 MeV) as 
given by Eqs. (14) and (15) for the case of i^C . The num- 
bers in parentheses correspond to the angular region with 
cos 6nv < —0.8. 





Nnp 


AT-An— i'nn 


^ ^ np 


jvrAnp— »nnp 


OPE 


1.00 (0.32) 


0.02 (0.00) 


0.91 (0.29) 


0.08 (0.02) 


OMEa 


0.89 (0.27) 


0.06 (0.02) 


0.76 (0.24) 


0.08 (0.02) 


OMEf 


0.87 (0.27) 


0.07 (0.02) 


0.7:! (0.2:!) 


0.08 (0.02) 



TABLE IV. Predictions for the ratio R-z = N„„/N„p for 
iHe and i^C. An energy thresholds of 30 MeV and two 
pair opening angle regions have been considered. The (pre- 
liminary) data are from KEK-E462 [27]. 





iHe 




i2ri 

A 






cos ^jvjv < —0.8 


all Onn 


cos Onn < —0.8 


all 6nn 


OPE 


0.25 


0.26 


0.24 


0.29 


OMEa 


0.51 


0.45 


0.39 


0.37 


OMEf 


0.61 


0.54 


0.43 


0.39 


EXP 


0.44 ±0.11 





TABLE V. Predictions for the weak interaction model in- 
dependent quantities NjPj^, N^j^ and N^j^ (integrated over 
all angles and for energies Tjv > 30 MeV) of Eqs. (14) and 
(15) for iHe and i^G. The numbers in parentheses correspond 
to the angular region with cos 9nn < —0.8. 





^IBn 


* nn 


jy2B 
' nn 


iHe 


0.84 (0.53) 


0.10 (0.02) 


0.54 (0.34) 




0.56 (0.30) 


0.27 (0.05) 


0.30 (0.12) 




* np 


^ np 


^ * np 


iHe 


0.20 (0.05) 


0.98 (0.49) 


0.55 (0.22) 




0.33 (0.08) 


1.22 (0.38) 


0.38 (0.11) 
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